Baryon kinetic energy loss in the color flux tube model. 

K.A. Lyakhov 1 

1 Frankfurt International Graduate School for Science, J.W. Goethe Universitat, 
Max-von-Laue Str. 1, D-60438 Frankfurt am Main, Germany 

Abstract 

This article generalizes Schwinger's mechanism for particles production in the arbitrary finite 
field volume. McLerran-Venugopolan(MV) model and iterative solution of DGLAP equation in the 
double leading log approximation for small x gluon distribution function were used to derive the 
new formula for initial chromofield energy density. This initial chromofield energy is distributed 
among color neutral clusters or strings of different length. This strings are stretched by receding 
nucleus. From the proposed mechanism of string fragmentation or color field decay based on exact 
solution of Dirac equation in the different finite volume, the new formulae for esimated baryon 
kinetic energy loss and rapidity spectrum of produced partons were derived. 

PACS numbers: 12.38.Mh, 25.75.Nq, 25.75.-q 



I. INTRODUCTION 



The process of partons production can be considered as the tunnelling across the ener- 
getic gap of width 2m± between the virtual energetic states inside Dirac sea to the real 
states-produced pairs with nonzero momentum. Semiclassical or WKB consideration of 
pair production was performed in [1], [18] and in many others. This article generalizes the 
Schwinger's mechanism for particles production in infinite field volume [3] for the case when 
caps of color flux tube recede from each other quasistatically. We have used low x behavior 
of parton distribution function to calculate dispersion of color charge per unit area and hence 
initial chromofield energy density. Moreover it will be shown how the total probability of 
string decay is related to the rapidity spectrum of produced partons. 

The current work uses the same methods as in [4], where finite field volume effects 
were considered first, and improves their results. Finite size effects on pair production in 
transverse direction were taken into account by applying MIT boundary conditions [? ]. The 
other possible but less general way to introduce influence of finite size effects on particles 
production is based on further development of Green functions method. For instance in 
[5] finite size effects were incorporated by expansion of Green functions on inverse volume 
occupied by field. 

At RHIC energies (7 ~ 100) nucleus can be represented as two massive sheets leaving the 
strong gluon field in their wake. It is instructive to split evolution of the system produced 
in heavy ion collisions on three characteristic stages. On the first stage immediately after 
collision the nuclei because of multiple soft gluons exchange acquire stochastic color charge. 
This color charges produce multitude of color flux tubes occupying space between receding 
streaks. Shortly afterwards flux tubes decay on prompt partons. In the second stage they 
lose part of their energy due to gluon radiation in cascade. And finally in the third stage 
secondary rescatterings drive the system to local thermal equilibrium. 

The subject of interest of this article will be the first stage of reaction where the screening 
of the color field (back reaction of produced plasma) may be disregarded. Parton distribution 
function calculated on this very initial stage can serve as initial condition for parton cascade 
model (PCM) [10], where particles are still energetic enough to apply pQCD for calculation 
of collision integrals and system can be treated as classical. 

In this paper probability for each vertex is calculated by solving wave equations where 
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FIG. 1: Energetic gap between positive and negative continuum in the presence of external field 
(linear string potential) as a function of coordinate. 

the volume occupied by field is restricted in transverse direction by the MIT boundary 
conditions and in the longitudinal direction by the distance L between colliding nuclei or 
capacitor plates [? ] . 



II. THE MODEL 



QCD Lagrangian and corresponding field equations read: 



£qcd{x) = ^!(x)[iYT) il -m}^!(x)-^F^ u F^ a (sum over repeated indices) a = 1 ... 8 (1) 



[i'fV VL -mYf!(x) = 0; 



V F a 



(2) 
(3) 

(4) 



dp-igs—Al. 

Gluon field is assumed gradually converting to plasma. This external field corresponds to 
the string potential = (A o (z),0). The energy stored in the color flux tube or gluon 

bag is regulated by its length L and its radius r , see below. Thus linear vector potential is 
given by 



A 



for z < (region I) 

— oz for < z < L (region II) 

—oL for z > L (region III) 



(5) 
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which is shown on Fig.l. 

In MIT bag model the energy stored in color flux tube is constrained by its volume AL, 
where A = nr^ is the string cross section. Therefore energy per unit length or string tension 
a can be expressed though the profile of the initial chromofield energy density €f(b, s) which 
will be estimated below: 

a = (e f (b,s) + B)A, (6) 

where B = Aq CD is the bag constant, and Aq C d = 200MeV. As it is shown on Fig. 2 the 
larger color flux tube radius the larger density of probability of pair production. This is 
easily explained by increasing string tension and therefore density of probability of string 
decay. 

Let us calculate initial chromofield energy density by following the ideas formulated in 
Refs. [2, 9] and widely used now as the Color Glass Condensate initial state. Within 
this picture, random color charges are generated on the nuclear sheets as a result of soft 
gluon exchange at the interpenetration stage. In a single event these charges fluctuate from 
point to point in the transverse plane. The charges also fluctuate from event to event, 
so that in average the areal charge is zero. It is convinient to introduce the color charge 
density, p(x~,s), as a function of coordinate in transverse plane and light cone variable 
x~ = (t — z)/y/{2). Following Refs. [2, 9] we assume it as a stochastic variable distributed 
with Gaussian weight: 



P[p] ~ exp 



- / d vd 2 s—^ — rTrp 2 (r],s) 
J 2/4 (?7, s) 



(7) 



where we have applied transformation of color charge density from light cone coordinate to 
pseudorapidity rj = \n(l/x~). In contrast to CGC model dependence on x~ is not important 
in MV model, and x~ dependence is reasonable to integrate out. Therefore areal charge 
density 

p(s) = J dr)p(r),s) (8) 
is random variable with gausssian distribution 

P[p] ~ exp (~ [ d 2 s] , (9) 
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Variance of color charge fluctuations in the nucleus a is denoted by p 2 a . These fluctuations 
are characterized by a certain scale in the transverse plane, which is related to the satura- 
tion scale a ^- introduced in high-density QCD [14]. For central Au+Au collisions at 



RHIC energies Q 2 xs lAAGeV 2 [15]. Thus, the "spots" on the transverse plane, where the 
color charge is essentially nonzero have the characteristic area s± = ira 2 = ^ ~ 0.09/m 2 . 
Since the transverse size of the baryonic slabs is much larger (about <tnn) this means that 
many string-like configurations(flux tubes) are stretched between the receding slabs. Each 
such configuration connects the spot of opposite charges, Qi, like in capacitor. In Abelian 
approximation the chromoelectric field strength in a flux tube is obtained from the Gauss 
theorem [8], 

E i = ^- = Pi (10) 

s± 

Then the force acting between opposite spots is 

1 1 

Fi = -QiEi = -p 2 s± = s±€i, (11) 

where e« is the energy density of the chromofield. It is important to note, that all flux tubes 
produce attractive force between opposite spots on the projectile and target slabs. Now let 
us divide the slabs into n small elements with equal transverse area s± so that they cover 
the total slab area a = ns±. The total force acting on each slab is then 



1 " 

= 9*l£p? (12) 



F 

2 ~ 

3=1 

Accordingly, we can represent the integral in eq.7 as the sum over all elements, so that 

P(PU ■ ■ ■ ,Pn) ~ exp (- 8 f ± I) = n exp (~ S -0) (13) 

\ i=l * / i=l ^ * ' 

Obviously the distribution of charge in each element follows the Gaussian distribution 



Therefore the mean force acting between two spots, as it follows from eqn.ll, is 

<Fi>=Y (15) 
The next step is to calculate the distribution of the total force acting between the slabs 
for an ensemble of events. This distribution is obtained by integrating 5(F — Fj) over all 

i 

charge densities pi with weight P(pi). for simplified case p\ = . . . = p? n = p 2 we have: 



w 



(f) = n (J p(Pi)d P )j s (f - l^izp^j 



V27r/x 2 / J V 2 / V /"V r ( n / 2 ) r(n/2)/x 2 \ /i 2 

(16) 

In the second expression we have used 0(4) symmetry of the inegrand and made transfor- 
mation to spherical coordinates in n-dimensional p space. As the result we get a gamma- 
distribution which has the following first moments: 

< f >= 

n 2 



^ = y-^ 2 (17) 

We see that the parameter ji 2 introduced in refs. [2, 9]in fact determines the mean force 
between the slabs and its dispersion. It is more convinient to express the mean energy 
density of the chromofield between the slabs: 

Taking into account that hard gluons and valence quarks act as source of chromofield we 
can calculate mean squared deviation of color charge per unit area of nucleus a as 

/i 2 = /i 2 , + /i 2 . (19) 

Quark color charge squared in a tube of transverse area d 2 s is the color charge squared per 
unit quark g 2 CV times the number of quarks in the tube dn q = N c N a (s, b)d 2 s 

= 9 2 C F ^ (20) 

where N a (s, b) is the number of participants from nucleus a at the given radius vector s 
and impact parameter b. Gluon color charge squared in a tube of transverse area d 2 s is the 
color charge squared per gluon g 2 N c times the number of gluons in the tube: 

i 



dn g = <N a (s,b) J G(x,Q 2 s )dx 

x 



> d 2 s, (21) 



where lower integration limit is defined from relation between Bjorken variable x and cm. 
beam energy ^Js: 

xo - ^, (22) 
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where Qq is the minimal transferee! four momentum. 

Low x gluon distribution function G(x, Q 2 ) is governed by the DGLAP equation [12]. In 
so-called double leading log approximation, where only terms proportional to InHnQ 2 are 
taken, it has the form 



G(x, Q 2 ) = —exp 

x 




U ~W c N f 



(23) 



where the value Ql = lGeV 2 was taken as a starting scale for Q 2 evolution. 
Running coupling constant is defined as g s (Q 2 ) = \f^a s (Q 2 ), where 

is the fine structure constant. 

Saturation scale inside nucleus a is obtained by the iterative solution of the following 
equation 

i 

Ql = ^\* s (Q 2 s )N a (b, s) J G(x, Q 2 )dx (25) 

xo 

Therefore in dense regime number of gluons per unit area is 

dn g _ C F Q 2 
d 2 s 4ir 2 a s 

Finally we have the following expression for mean color charge squared: 



(26) 



/i 2 = C F N C [ Ana s N a + ^ ) (27) 



III. PROBABILITY OF STRING DECAY 

In this section we calculate total probability of string decay. Complicated further evo- 
lution of cascade related with numerous branchings of secondaries is out of scope of our 
consideration. Due to the three possible interactions which follow from C-qcd string can 
decay into quark-antiquark, gluon pair or tree gluons. Then the total density of probability 
is 

W tot (p, z, L) = W q g(p, z, L) + W gg (p, z, L) + W ggg (p, z, L) (28) 
Last term will be disregarded due to lack of analytic solution for 3-body problem. 
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The density of probability of gluon pair production is calculated by the formula 



Wggip, Z, L) = VgWg(p, Z, L) , 



(29) 



where v g = 2(N 2 — 1); W g — ^*^, and ^ is solution of Klein-Gordon equation for massless 
particles (55). 

The density of probability of quark- ant iquark pair production is defined as 



W q g(p, Z, L) = VqWqq{p, Z, L) , 

where v q = N c Nf; W qq - = and * is solution of Dirac equation (34). 
Solution of the Dirac equation can be represented by the following series: 

V(x) = * e MKPxX + VyV ~ Et)], 
V 

where 

r=l,2 

Eigenvectors ty r are Dirac spinors corresponding to the different spins [20]: 

( ^1 + \ 

A r ^ 2 - i>i 

\ -B r (A r fa-fa)/ 
making of the eigenfunctions of the squared Dirac equation: 



(30) 



(31) 



(32) 



\l/ r = a r 



(33) 
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i— ) — ill 2 ::L l(T I ; ',. = 0. 



(34) 



Cylindrical boundary conditions applied to the bag surface discretize transverse momentum 
inside the bag. The ground state of transverse momentum p\ — ci/r is included to effective 
parton mass m t [? ]; m t — >/m 2 + (p° ± ) 2 , where m is the current parton mass; c\ = 1.4347. 
The spin coefficients were calculated in [20]: 



A, 



mj_ ± ip± 



±i 



.Px + iPy 



(35) 



m± p± 

where m± = \fm 2 + p 2 ± ; p± = ^Jp 2 + p 2 . Therefore the decomposition (32) in accordance 
with (33) is transformed to: 

§ = /ii^i + /i 2 ^2, (36) 
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where [i\ = 2cti 



/l\ 

1 

w 



; fJ>2 = 2a 2 



/ J2H \ 

Vx+ipy 

mj_ 

m t 
mj_ 

Px+iPy 
mi 



Spinor normalization constants are defined from condition /4/i r = 1: a r = r = 1,2. 
After separation of variables (34) takes Schrodinger-type form: 



where 



Pr 



p\ for z < (region I) 

(oj{p) + crz) 2 — m\ ±icr for < z < L (region II) 
p\ for z > L (region III), 



(37) 



(38) 



where E L = u){p); E R = uj{p) + <jL; p R;L = y E R L — m\ and u(p) < —m±. Solutions in 
regions I, III are the following: 



Ie ik,r + R ( L yk R r for z < o ( region I) 

T(L)e ikTT for z > L (region III), 



(39) 



where R, T are amplitudes of the reflected, and transmitted waves, respectively and four- 
momenta in this regions are 

fc/ = (E L ,p x ,p y ,-p L ) (40) 

k R = (E L ,p x ,p y ,p L ) (41) 

k T = {E R ,p x ,p y ,p R ) (42) 

kr = E — p x x — p y y - p z z (43) 
Solution in the region II is represented as a sum of the linear independent parabolic cylinder 

2 

function of order v r = ip r — 1/2; p r = — ^ ± s; s = | [21]: 

^ r = A(L)£>° r (e^CM) + B(L)Dl(e^a P ,z)) (44) 

where CO ?- 2 ) — \J~^( uj (p) + <T - 2 )- Coefficients A(L),B(L) are found by matching the wave 
function at the boundaries. The boundary conditions implement the wave function conti- 
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nuity at z = and z = L: 

' (I + R(L))jj r = (A(L)D° L + B(L)Dl)^ 
ip L (-I + R(L))fi r = {A{L)D« + B(L)D^)fi r 
T(L)e ipRL fi r = (A(L)D° R + B(L)D R )fi r 
( i PR T(L)e^ L fi r = (A(L)D° R + B(L)D R )fj, r , 
where D ^ = Dl^(E LjR e tJ z), and amplitude of incident wave is / = 1. Therefore 

2 D% - ip L D R 



(45) 



B(L) 



qi i PR (D R -qD R ) + qD R '-D 

A(L) = -- B(L)q; 
Qi 



R 



R{L) = B{L){D l L -qDl) + -Dl-l- 



T(L) = e~ ipRL 



where 



B(L)(D R -qD R ) + -D R 



D o> D}-^- 



q 1 = D° L --±;q = 

IPL 



Gluon production is described by Klein-Gordon equation: 
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or 



Solution of this equation is the following: 

* = Yl Yl * exp[i(PxX + VyV - Et)]e x , 

V A=±l 

where t\ polarizations of vector particles. 

After separating of variables (51) takes the Schrodinger-type form: 

d 2 



dz 2 



p A ^ = 0, 



where 



V 



(46) 
(47) 
(48) 
(49) 

(50) 
(51) 
(52) 



p\ for z < (region I) 

(u(p) + az) 2 - m 2 ± for < z < L (region II) 
p\ for z > L (region III), 

and m_|_ = \J{f^) 2 + v\ is the gluon transverse mass. In the region II solution 



equation is 



where p = — ^ 



* = A(L)D»(e^((p, z)) + B(L)Dl(e^C(p, z)) 



(53) 



(54) 



of this 



(55) 
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IV. TRANSMISSION AMPLITUDE OF PAIR PRODUCTION IN INFINITE 
VOLUME 



In this section we will show that in infinite volume occupied by field the density of 
probability or transmission amplitude has the proper Schwinger asymptotic. Transmission 
amplitude is expressed by the formula (49) which is transformed to the following 

2p L (D° R D 1 R -D° R D 1 R )e- ipRL 

At large arguments parabolic cylinder function [21] takes the following form D R — > e ^~ , 
D R \ ±ip R ^ L D R 1 L . Therefore asymptotic of transmission amplitude reads: 

^SHi— - 

From latter it follows that amplitude of pair production in infinite field volume approaches 
the following limit: 

w = lim |T(L)| 2 ~ exp ( . (58) 
L^oo y cr J 

V. PARTONS RAPIDITY DISTRIBUTION AND EXPANDING CHRO- 
MOFIELDS 

In the initial stage of reaction all the kinetic energy lost by baryons is transformed to 
chromofield: 

E f {b) = J d 2 s(e f (b, s) + B)z(b, s), (59) 
where average string length is 

oo oo 
oo „ „ 

*(M) = J2 d 'P dzzW tot (p,z,L), (60) 

T — T ■ J 

where L min = ^ =L . On the other hand finally all chromofield energy is transformed to 
partonic plasma 

E f (b) = J d 3 P T(p) = J d 3 pu(p)^- (61) 

Therefore rapidity distribution of partons generated in the slab-slab collision is calculated 
by the formula: 

oo 

^-(b) = J d 2 p ± r( P ) (62) 
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FIG. 2: Density of probability of quark-antiquark pair production Wgg(z,p±, y, L) as a function of 
coordinate inside the color flux tube of length L = 10/m at the fixed rapidity y = 0.21 and average 
transverse momentum < p± >= 0.5GeV. 

VI. CONCLUSIONS 

The results of this article are following: l)the total density of probability of partons 
production is obtained by the exact solutions of squared Dirac and Klein-Gordon equations in 
the linear vector potential; 2) Rapidity spectra of the field energy and produced particles are 
calculated from the energy conservation law; 3)It was shown that the density of probability 
of particles production in the infinite field volume approaches the classical Schwinger's result. 

Numerical simulations were performed for most central Au+Au collisions with cm. energy 
y/s = 200AGeV . The coordinate dependencies of a gg-pair production density of probability 
at the different color flux tubes radii (r = 0.25 fm, 0.5 fm, lfm) are shown on Fig.2. Quarks 
have current mass m q ~ 8MeV and gluons are assumed massless. 
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